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On  the  basis  of  general  relationships  of  the  thin  body  theory,  the  paper  /3* 
investigates  the  flow  past  a thin  plane  wing  of  low  aspect  ratio  in  a wind 
tunnel  with  a circular  test  section.  The  corrections  to  the  lift  coefficients 
and  position  of  the  wing  pressure  center  are  determined  by  allowing  for  the 
boundary  constraint  effect.  The  solutions  obtained  are  compared  with  known 
formulas  for  experimental  corrections. 

We  will  consider  the  motion  of  a thin  plane  wing  of  low  aspect  ratio  in 
a cylindrical  wind  tunnel  of  circular  cross  section  (Fig.  1).  In  accordance 


with  the  thin  body  theory,  the  three-dimensional  flow  in  a wind  tunnel  caused 
by  the  motion  of  the  wing  will  be  replaced  by  a two-dimensional  flow  in  a 
transverse  plane  fixed  in  the  tunnel.  The  flow  in  this  plane  is  caused  by  ver- 
tical displacement  of  the  plate  at  velocity  Vq  = V^a  (wing  cross  section  wake) 


and  by  a change  in  its  width.  This  flow  will  result  in  a change  in  the  apparent 
mass  of  the  plate,  proportional  to  the  increment  of  the  wing  lift  on  the  length 

dx,  = V dt: 

1 «> 


dY  _ _ ,,2  dm 

dx,  dt  ~~  ~a  dx, 


(J) 


If  (1)  is  summed  up  over  the  length  of  the  root  chord  (0  ? x,  ^ 1),  the  total 
wing  lift  is  obtained:  1 

y — Vl'O.m  (1),  (2) 

where  m(l)  is  the  value  of  the  apparent  mass  in  section  x,  •*  1 , coinciding  with 
the  maximum  wingspan. 
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Using  (1),  one  can  determine  the  longitudinal  wing  moment  relative  to  the 
o^z^  axis  passing  through  the  wing  apex: 


C dY 
- J dxxl 


We  will  introduce  the  lift  and  moment  coefficient 


; mt 


and  determine  the  position  of  the  wing  pressure  center  relative  to  the  apex  in 

mz 

fractions  of  the  root  chord  as  x = - — . 

P Cy 

• i 

From  (2)  and  (3),  we  obtain  the  following  formulas: 

cVa,2ir:  xp=  ,“S7n{mwrfx-  (4) 

According  to  (4),  to  determine  the  lift  coefficient  and  position  of  the  pressure 
center,  it  is  necessary  to  know  the  chordwise  distribution  of  the  apparent  mass 
of  the  wing  cross  sections. 

To  calculate  the  apparent  mass  of  the  plate  moving  perpendicularly  to  its 
plane  in  the  region  bounded  by  the  circle,  we  will  distribute  a vortex  sheet 
of  variable  strength  y(z)dz  over  the  length  of  the  plate.  A vortex  element  of 
y(z)dz  will  be  placed  at  an  arbitrary  point  of  the  plate  -SL  4 z 4 1.  This  vortex 


produces  a vertical  velocity  of  -cL  yiz)dz  at  a fixed  point  of  the  plate  S.  T1 

action  of  all  the  vortices  is  determined  by  summing  up  over  the  length  of  the 
plate 

J_  f i (?)  dz 

2"  J * “ * • 

-/ 

The  influence  of  the  flow  boundaries  will  be  taken  into  account  by  means 
of  an  additional  Vertex  with  a strength  of  ±y(z)dz  (the  plus  sign  refers  to 
a wind  tunnel  with  an  open  test  section,  and  the  minus  sign,  to  a closed  test 

r2 

section),  placed  at  a distance  z^  = — from  the  tunnel  center  (Fig.  1). 

The  effect  of  the  wind  tunnel  on  the  plate  is  made  up  of  the  effects  of 
additional  vortices  and  is  expressed  by  the  fact  that  at  point  s of 
the  wing,  an  additional  vertical  velocity  is  obtained,  equal  to 


at  a fixed  point  of  the  plate  S.  The 


_1_  p t (?)  dz 
x 2k  \ W 

iT-  . 


Since  the  total  velocity  at  point  s should  be  equal  to  the  displacement 
velocity  of  the  plate  VQ,  then 

(5) 


I 


Here  k = ±1,  where  the  plus  sign  refers  to  a tunnel  with  an  open  test  section 
(free  jet),  and  the  minus  sign,  to  a tunnel  with  a closed  test  section  (solid 
walls). 

Relation  (5)  is  an  integral  equation  for  determining  the  unknown  strength 

Y(z)» 

We  will  switch  to  new  variables  in  Eq.  (5)  after  performing  the  substitution 
z *>  £cos  $ ; s = £cos6;  y(z)  y($);  0 < $ it.  After  some  simple  transformations, 

we  obtain  (p  = ^) 


1C  * 

^ (»)  sin  , . f t(!>)  cos»sin»d»_  « w 

\ cos  » — cos  6 ' \ 1 „ 0 

J J -,-cos»cos0 


Equation  (5)  or  (6)  has  an  exact  solution  corresponding  to  the  motion  of  the 
wing  in  an  unbounded  flow.  In  this  case  R = °°,  and  the  second  term  in  Eqs.  (5) 
and  (6),  which  allows  for  the  boundary  effect,  vanishes,  and  Eq.  (5)  becomes 


Inverting  this  integral  on  the  condition  of  unboundedness  of  y (z)  on  the  ends 
of  the  [-£,  £]  interval,  we  obtain 

T-  “ y^rz?  = 2V.°tan  ^ (7) 

To  find  the  solution  of  Eqs.  (6)  for  finite  wind  tunnel  radii  R,  the  unknown  / 

function  y($)  will  be  .represented  in  the  form  of  the  following  series  with  un- 
known coefficients: 

7 (0)  = 2Vt  j^o ctg 0 + 2 Ain  sin  n0 j . . (8) 

Here  the  first  term  for  Aq  = 1 is  the  solution  in  an  unbounded  flow,  and  the 
second  term,  when  A^  4 1,  determines  the  boundary  effect. 

If  the  coefficients  A^  (n  = 0,  1*  2,  ...)  are  known,  the  apparent  mass  of 
the  plate  is  determined  according  to  Ref.  1 by 

m = sP/*(/l,  + £).  (9) 

In  the  special  case  of  R = “ (unbounded  flow),  AQ  - 1,  = 0,  and 

« it (10) 

lo  determine  A^  (n  =0,  1,  ...)  in  the  general  case  with  R ^ °°,  we  will  substitute 
(8)  into  (6).  After  a caries  of  transformations,  we  obtain 


I H 


^ v / . , ' . / r — : Ain  fees  2/i0  ■ 

» I — |»'ct)0;3  vl  + j I— n-cos-«/j  a~| 
h yXM-OcwX*— >'0  1 
(i  y'i  — (i*  cow  w j 


It  a scries  of  values  0 = -^s  (s  = 0,  1 p)  is  fixed  in  (11),  a system 

r>t  equations  is  obtained  for  determining  the  coefficients  of  the  series  A^, 

,,,,  • The  remaining  coefficients  must  be  set  equal  to  zero. 

This  method  was  used  to  check  calculations  of  the  coefficients  for  different 
numbers  of  fixed  points,  the  analysis  of  which  showed  that  the  systems  of  alge- 
braic equations  are  very  stable,  and  the  apparent  mass  values  calculated  by 
kiitterent  .approximations  already  practically  coincide  starting  with  p £ 2. 


MB 


FIGURE  2 

The  results  obtained  from  calculations  of  the  relative  magnitude  of  the  /7 
apparent  mass  of  the  plate  (m*  = — ) for  p = 6 are  given  in  Fig.  2 (dashed  curve) . 

IDoo 

If  we  confine  ourselves  to  the  condition  y«  1,  we  can  obtain  an  approxi- 
mate analytical  solution  of  Eq.  (5)  or  (6). 

The  denominator  of  the  second  term  of  Eq.  (5)  will  be  represented  in  the  form 
! = i ! ~_i/i  i « i’****  , ''  \ » 

£ K* 

* * — 

Tills  being  taken  into  account,  Eq.(5)  is  written  as 

i i 

f t (r)  dr  , k C , . , 

\-TT7  + *5  U <*)**•  (12) 

It  we  now  turn  to  the  new  variable  $ (z  = Hcos^)  and  use  expansion  (8),  we  ob- 
tain from  (12)  in  the  accepted  formulation 

^ — I Ain  e=0(rt  «=  1,  2,  . . .).  (13) 

cording  to  (9)*  the  relative  apparent  mass  will  be  written  in  this  approximation 


Comparison  of  calculations  made  by  means  of  formula  (14)  with  the  exact  solu- 
tion (Fig.  2)  shows  their  satisfactory  agreement  for  y ^ 0.8. 

The  approximate  solution  (14)  is  used  below  for  analyzing  the  effect  of 
the  tunnel  walls  on  the  aerodynamic  characteristics  of  a low-aspect-ratio  wing. 

The  effect  of  the  tunnel  walls  on  the  lift  coefficient  is  established 

directly  from  (4)  and  (14).  We  will  denote  by  and  the  values  of  the 

*00  *£ 

derivatives  of  the  wing  in  an  unbounded  flow  and  in  a wind  tunnel.  The  rela- 
tionship between  them  will  be  as  follows: 

c*  = — c*  =cf.c* 

m*  't  t Yt 


CL  = (!  +1  v-1)0*? 


It  is  interesting  to  note  that  the  correction  factor  in  formula  (15)  is  in 
complete  agreement  with  the  existing  corrections  in  a circular  tunnel2  if  in 
the  latter  corrections  one  takes  A -*•  0. 

/8 

To  allow  for  the  effect  of  the  tunnel  walls  on  the  position  of  the  wing  pressure- 
center,  it  is  necessary,  according  to  (4),  to  know  the  distribution  of  the 
apparent  mass  over  the  length  of  the  wing.  This  distribution  depends  on  the 
planform  of  the  wing. 

We  will  consider  a class  of  wings  whose  planform  is  given  by  the  equation 

(*i)  = l (*i)  = Lx»  (0  < * < I), 


where  I.  is  the  wing  semispan. 

The  change  in  apparent  mass  along  the  length  of  the  wing  is  determined 
from  (14): 

m (a)  «p  D . (17) 

1 ±2 

The  position  of  the  pressure  center  of  the  wing  relative  to  its  apex  in 
fractions  of  the  root  chord  is  calculated  from  formula  (4): 


-0+H*  J-r 

1 T o 


+ 2 M5*2p 


In  an  unbounded  flow  y = 0,  and  from  (18)  we  obtain 


'p-  = 

V 

X 

uv  Ini  induce  rhc  relative  quantity  6 - P—  . We  then  obtain  the  following  rela- 

Xpa. 

tlon  between  the  values  of  the  pressure  center  in  an  unbounded  flow  x and  in 

p» 

® wind  luiniel  x : 

P. 


where  according  to  (17)  and  (18),  6 is  determined  as  follows: 

2 p 


5 = 


2 P+  I 


. v 'J  i + j^p 


(21) 


The  magnitude  of  the  correction  factor  6 for  certain  planforms  (p  values)  of 
wings  is  shown  in  Fig.  3. 
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FIGURE  3 


The  results  of  the  calculations  show  that  the  effect  of  the  tunnel  walls 
on  the  position  of  the  pressure  center  is  less  pronounced  than  that  on  the 

magnitude  of  the  derivative  C^- 
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On  the  basis  of  general  relationships  of  the  thin  body  theory,  the  paper  /3* 
investigates  the  flow  past  a thin  plane  wing  of  low  aspect  ratio  in  a wind 
tunnel  with  a circular  test  section.  The  corrections  to  the  lift  coefficients 
and  position  of  the  wing  pressure  center  are  determined  by  allowing  for  the 
boundary  constraint  effect.  The  solutions  obtained  are  compared  with  known 
formulas  for  experimental  corrections. 

We  will  consider  the  motion  of  a thin  plane  wing  of  low  aspect  ratio  in 
a cylindrical  wind  tunnel  of  circular  cross  section  (Fig.  1).  In  accordance 


FIGURE  1. 


with  the  thin  body  theory,  the  three-dimensional  flow  in  a wind  tunnel  caused 
by  the  motion  of  the  wing  will  be  replaced  by  a two-dimensional  flow  in  a 
transverse  plane  fixed  in  the  tunnel.  The  flow  in  this  plane  is  caused  by  ver- 
tical displacement  of  the  plate  at  velocity  Vn  - V a (wing  cross  section  wake) 

u 00 


and  by  a change  in  its  width.  This  flow  will  result  in  a change  in  the  apparent 

mass  of  the  plate,  proportional  to  the  increment  of  the  wing  lift  on  the  length 

dx.  = V dt: 

1 <“ 


dV  _ d y2  dm 

dx,  dt  ~ *>a  dx,  ‘ 


(1) 


If  (])  is  summed  up  over  the  length  of  the  root  chord  (0  £ x.  $ 1),  the  total 
wing  lift  is  obtained:  1 

Y = VU.m.{  1),  (2) 

where  rn(l)  is  the  value  of  the  apparent  mass  in  section  x^  *»  1,  coinciding  with 
the  maximum  wingspan. 


•k 
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Using  (1) , one  can  determine  the  longitudinal 
axis  passing  through  the  wing  apex: 

i 

M*>  1=1  j Xi  !xldx‘ 

0 


wing  moment  relative  to  the 


<3) 


We  will  introduce  the  lift  and  moment  coefficient 


and  determine  the  position  of  the  wing  pressure  center  relative  to  the  apex  in 

mz 

fractions  of  the  root  chord  as  x = - — . 

P cY 

From  (2)  and  (3),  we  obtain  the  following  formulas:^- 


/A 


c%  = 27T- : *p=  1 — m Mdx-  <4) 

According  to  (A),  to  determine  the  lift  coefficient  and  position  of  the  pressure 
center,  it  is  necessary  to  know  the  chordwise  distribution  of  the  apparent  mass 
of  the  wing  cross  sections. 

To  calculate  the  apparent  mass  of  the  plate  moving  perpendicularly  to  its 
plane  in  the  region  bounded  by  the  circle,  we  will  distribute  a vortex  sheet 
of  variable  strength  y(z)dz  over  the  length  of  the  plate.  A vortex  element  of 
y(z)dz  will  be  placed  at  an  arbitrary  point  of  the  plate  -l  4 z 4 l.  This  vortex 

« 1 y(z)dz 

produces  a vertical  velocity  of  — — - at  a fixed  point  of  the  plate  S.  The 

action  of  all  the  vortices  is  determined  by  summing  up  over  the  length  of  the 
plate 

2it 


z — $ ■ 


-/ 


The  influence  of  the  flow  boundaries  will  be  taken  into  account  by  means 
of  an  additional  vortex  with  a strength  of  ±y(z)dz  (the  plus  sign  refers  to 
a wind  tunnel  with  an  open  test  section,  and  the  minus  sign,  to  a closed  test 

r2 

section),  placed  at  a distance  z^  = — from  the  tunnel  center  (Fig.  1). 

The  effect  of  the  wind  tunnel  on  the  plate  is  made  up  of  the  effects  of 
all  the  additional  vortices  and  is  expressed  by  the  fact  that  at  point  s of 
the  wing,  an  additional  vertical  velocity  is  obtained,  equal  to 


i f Liilif 
2,  \ £ • 


-i 


Since  the  total  velocity  at  point  s should  be  equal  to  the  displacement 


veJocity  of  the  plate  VQ,  then 


f£r  + * fM 

-i  i,  > 


1 (')  di 


(5) 


/5 


— s 


Here  k = ±1,  where  the  plus  sign  refers  to  a tunnel  with  an  open  test  section 
(free  jet),  and  the  minus  sign,  to  a tunnel  with  a closed  test  section  (solid 
walls). 

Relation  (5)  is  an  integral  equation  for  determining  the  unknown  strength 

Y(z). 

We  will  switch  to  new  variables  in  Eq.  (5)  after  performing  the  substitution 
z = £cos  $ ; s = £cos0;  y(z)  y($);  0 ^ ^ it.  After  some  simple  transformations, 

£ 

we  obtain  (y  = — ) 


(*  T (»)  sin  . f y (»)  cos  » sin  „„ 

\ cos  i>  — cos  6 ' \ 1 „ 0 

J . J ~j-cos»cos0 


Equation  (5)  or  (6)  has  an  exact  solution  corresponding  to  the  motion  of  the 
wing  in  an  unbounded  flow.  In  this  case  R = °°,  and  the  second  term  in  Eqs.  (5) 
and  (6),  which  allows  for  the  boundary  effect,  vanishes,  and  Eq.  (5)  becomes 


j*  7-  Wdz  _ 


Inverting  this  integral  on  the  condition  of  unboundedness  of  y (z)  on  the  ends 
of  the  [-£,  £]  interval,  we  obtain 

7-  = = 2v'ocotan  S'  (7) 

To  find  the  solution  of  Eqs.  (6)  for  finite  wind  tunnel  radii  R,  the  unknown  j 

function  y($)  will  be  represented  in  the  form  of  the  following  series  with  un- 
known coefficients: 

T (0)  = 2V0  J\z40ctgO  + 2 Ain  sin  n°]  • 

Here  the  first  term  for  Aq  = 1 is  the  solution  in  an  unbounded  flow,  and  the 
second  term,  when  A^  4-  1,  determines  the  boundary  effect. 

If  the  coefficients  A^  (n  = 0,  1,  2,  ...)  are  known,  the  apparent  mass  of 
the  plate  is  determined  according  to  Ref.  1 by 

m = + (9) 

In  the  special  case  of  R = » (unbounded  flow),  Aq  = 1,  - 0,  and 

m..  = rp /*.  (10) 

To  determine  A^  (n  = 0,  1,  ...)  in  the  general  case  with  R 4 we  will  substitute 
(8)  "In^c  (6).  After  a series  of  transformations,  we  obtain 


[ 1 *1*  " 2 Aln  [c°s2«a  — 

) 1 — i*1  iu»#  0(1  + ) I — /a4  rus*  0)  J n~  I 

(l  -f  V I — I*4  COS*  'I  ) ,n\ 


• • • t 


If  a series  of  values  0 = -^pS  (s  = 0,  1,  . . . , p)  is  fixed  in  (11),  a system 
of  equations  is  obtained  for  determining  the  coefficients  of  the  series  A^, 

A, 


A . The  remaining  coefficients  must  be  set  equal  to  zero. 

2p 


Tills  method  was  used  to  check  calculations  of  the  coefficients  for  different 
numbers  of  fixed  points,  the  analysis  of  which  showed  that  the  systems  of  alge- 


braic equations  are  very  stable,  and  the  apparent  mass  values  calculated  by 
dltterent  approximations  already  practically  coincide  starting  with  p ^ 2. 
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The  results  obtained  from  calculations  of  the  relative  magnitude  of  the 
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apparent  mass  of  the  plate  (m.  = — ) for  p = 6 are  given  in  Fig.  2 (dashed  curve). 
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If  we  confine  ourselves  to  the  condition  y«  1,  we  can  obtain  an  approxi- 
mate analytical  solution  of  Eq.  (5)  or  (6). 


The  denominator  of  the  second  term  of  Eq.  (5)  will  be  represented  in  the  form 

I 2 I 


A" 


7?r+ 


This  being  taken  into  account,  Eq. (5)  is  written  as 


P t (?)  ch  k p . , 

jT=T+/fi 


(12) 


If  we  now  turn  to  the  new  variable  (z  =-  £cos$-)  and  use  expansion  (8),  we  ob- 
tain from  (12)  in  the  accepted  formulation 


/10  = 


l + Tt*’ 


•I  A = 0 (n  « !,  2.  . . .). 


02} 


According  to  (9),  the  relative  apparent  mass  will  be  written  in  this  approximation 
as 

I 


Ml 


(H) 


H 2 il* 


S6 


Comparison  of  calculations  made  by  means  of  formula  (14)  with  the  exact  solu- 
tion (Fig.  2)  shows  their  satisfactory  agreement  for  y £ 0.8. 

The  approximate  solution  (14)  is  used  below  for  analyzing  the  effect  of 
the  tunnel  walls  on  the  aerodynamic  characteristics  of  a low-aspect-ratio  wing. 

The  effect  of  the  tunnel  walls  on  the  lift  coefficient  is  established 

Ct  ct 

directly  from  (4)  and  (14) . We  will  denote  by  C and  C the  values  of  the 

Y°°  Yt 

derivatives  of  the  wing  in  an  unbounded  flow  and  in  a wind  tunnel.  The  rela- 
tionship between  them  will  be  as  follows: 

1 


or 


- nr,  Cvt  ~ Ct crt‘ 
CL  = (*  ^ 2.  ^ ) Cyf 


(15) 


It  is  interesting  to  note  that  the  correction  factor  in  formula  (15)  is  in 
complete  agreement  with  the  existing  corrections  in  a circular  tunnel2  if  in 
the  latter  corrections  one  takes  A 0. 

/ 8 

To  allow  for  the  effect  of  the  tunnel  walls  on  the  position  of  the  wing  pressure 
center,  it  is  necessary,  according  to  (4),  to  know  the  distribution  of  the 
apparent  mass  over  the  length  of  the  wing.  This  distribution  depends  on  the 
planform  of  the  wing. 

We  will  consider  a class  of  wings  whose  planform  is  given  by  the  equation 

(*i)  = l (Xj)  =zLx*  (0  < * < 1), 


(16) 


where  I.  is  the  wing  semispan. 


The  change  in  apparent  mass  along  the  length  of  the  wing  is  determined 


from  (14): 


m (a)  = :tp L2 


r IP 


» ±y  -i »*2p 


(17) 


The  position  of  the  pressure  center  of  the  wing  relative  to  its  apex  in 
fractions  of  the  root  chord  is  calculated  from  formula  (4): 


(16) 


In  an  unbounded  flow  y = 0,  and  from  (18)  v?e  obtain 

2P 


Vp+l- 


(19) 


*»e  introduce  the  relative  quantity  6 - — . We  then  obtain  the  following  rela- 

XP™ 

f °0  • 

tion  between  the  values  of  the.  pressure  center  in  an  unbounded  flow  x and  in 

Poo 

a wind  tunnel  x : 


' 6 'pi 


't 


where  according  to  (17)  and  (18),  6 is  determined  as  follows: 


2 p 


The  magnitude  of  the  correction  factor  <5  for  certain  planforms  (p  values)  of 
wings  is  shown  in  Fig.  3. 


FIGURE  3 


The  results  of  the  calculations  show  that  the  effect  of  the  tunnel  walls 
on  the  position  of  the  pressure  center  is  less  pronounced  than  that  on  the 

magnitude  of  the  derivative 
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